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\ Abstract. Using the loop equations we find an explicit expression for genus 1 correction in 

I hermitian two-matrix model in terms of holomorphic objects associated to spectral curve arising 

J> ■ in large N limit. Our result generalises known expression for in hermitian one-matrix model. 

. We discuss the relationship between , Bergmann tau-function on Hurwitz spaces, G- function of 
Frobenius manifolds and determinant of Laplacian over spectral curve. 

o 

\ In this letter we derive an explicit formula for the l/N"^ correction to free energy F of hermitian 

I two-matrix model: 

! e~^^^ 

Oh 



It is hard to overestimate the interest to random matrix models in modern physics and mathematics; 
we just mention their appearance in statistical physics, condensed matter and 2d quantum gravity (see 
e.g. PP) . The expansion F = ^'q^qN~'^^F^ [N is the matrix size) in hermitian matrix models is 
^ \ one of the cornerstones of the theory, due to its clear physical interpretation as topological expansion 
of the functional integral, which appears in — > oo limit; in statistical physics the term F^ plays 
the role of free energy for statistical physics model on genus G Riemann surface. From the whole 
zoo if the random matrices one of the simplest is the hermitian one-matrix model with partition 
function ^ — J dMe~^^^'^^^'^'> (V is a polynomial), which can be used as testing ground for the 
methods applied in more general situations of two- and multi- matrix models. The most rigorous 
way to compute the 1/A^ expansion for both one-matrix and two-matrix models is based on the loop 
equations. The loop equations follow from the reparametrization invariance of matrix integrals; for 
one-matrix case the loop equations were used to compute F^ (see 0). Later the loop equations were 
written down for the case of two-matrix model |31 11] and F^ was found for the case when the spectral 
curve has genus zero and one j^; for arbitrary genus of spectral curve of two-matrix model only the 
leading term F^ is known (see |31). 

Let us write down polynomials Vi and V2 in the form Vi{x) = Yl!k=i it^^ ^^{y) = Z]fc=L^ X^^^- 
It is sometimes convenient to think of Vi and V2 as infinite formal power expansions: Vi{x) = 
X^fcLi x^*^' ^2(2/) = SfcLi ^y'^j where coefficients Ufc vanish for k > di + 2, and Vk vanish for 
k > d2 + 2. According to this point of view the operators of differentiation with respect to coefficients 
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of Vi and V2 have the fohowing meaning (see 0): 



6Vi{x) x' I £^ \ Mfc=o ,fc>di+2 ' <^V2(y) 



(2) 

As it was discussed in detail in 0, Q is a formal notation which makes sense only order by order in 
the infinite power series expansion; it allows to write an infinite number of equations at once. Consider 



the resolvents W(x) = x-M-^ / ^(y) = 'Ny''^ y-M2 / "^^^ ^^^^ energy of two-matrix model 
satisfies the following equations with respect to coefficients of polynomials Vi and V2: 

SF , , 6F ~ , 



6Vi{x) ' ' ' 6V2{y) 

The equations Q were solved in in the zeroth order assuming the finite-gap structure of distribution 
of eigenvalues of Mi (and, a posteriori, also of M2) as ^ 00. Here we find the next coefficient F^, 
using the loop equations. The spectral curve C is defined by the following equation, which arises in 
the zeroth order approximation: 

£'>{x, y) := (Fi'(x) - y){V2'{y) - x) - V\x, y) + 1 = (4) 

where the polynomial of two variables V^{x,y) is the zeroth order term in expansion of the 

polynomial 

the point P € C of the curve is the pair of complex numbers {x,y) satisfying @ (on the "physical" 
sheet the equation of spectral curve ^ defines an implicit function y{x), which gives the zeroth 
order approximation to V2'{x) — W{x)). The spectral curve ^ comes together with two meromorphic 
functions /(P) = x and g{P) = y, which project it down to x and y-planes, respectively. These 
functions have poles only at two points of £., called 00/ and oo^: at 00/ function f{P) has simple 
pole, and function g{P) - pole of order di with singular part equal to Vi'{f{P)). At oo^ the function 
g{P) has simple pole, and function f{P) - pole of order d2 with singular part equal to V2'{g{P)). 
In addition, in our normalization of partition function we have the asymptotics VV(x) ~^^oo 
1/x + . . . and VV(y) ~y^oo 1/y + • • • , which imply Kesoofgdf = 1 and Resoogfdg = 1, respectively. 
Therefore, one gets the moduli space Ai of triples {jC,f,g), where functions / and g have this pole 
structure. The natural coordinates on this moduli space are coefficients of polynomials Vi and V2 
and g numbers, called "filling fractions" ea = 2^ £j gdf, where {aa,ba) is some basis of canonical 
cycles on £. The additional constraints which should be imposed a posteriori to make the "filling 
fractions" dependent on coefficients of polynomials Vi and V2 are (according to one-matrix model 
experience, these conditions correspond to non-tunneling between different intervals of eigenvalues 
support): gdf = 0. Denote the zeros of differential d/ by Pi, ... , Pmi {rni = d2 + 2g + 1) (these 
points play the role of ramification points if we realize L as branched covering by function f(P)); 
their projections on x-plane are the branch points, which we denote we denote by Ij := f{Pj) ■ The 
zeros of the differential dg (the ramification points if we consider C as covering defined by function 
g{P)) we denote by Qi, . . . , Qm2 ("^-2 = di + 2g + 1); there projections on y-plane (the branch points) 
we denote by fij := g{Qj)- We shall assume hat our potentials Vi and V2 are generic i.e. all zeros 
of differentials df and dg are simple, and none of the zeros of df coincides with a zero of dg. If is 
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well-known (see for instance jS]) how to express all standard algebro-geometrical objects on L in terms 
of the previous data. In particular, the Bergmann bidifferential B{P,Q) = dpdQliiE{P,Q) {E{P,Q) 
is the prime-form) can be represented as follows: 

The Bergmann bidifferential has the following behaviour near diagonal P ^ Q: B{P,Q) = {{z{P) — 
z{Q))-^ + ISb{P) + o{l)}dz{P)dz{Q), where z{P) is some local coordinate; SsiP) is the Bergmann 
projective connection . Consider also the four-differential D{P, Q) = dpdg In E(P, Q), which has on the 
diagonal the pole of 4th degree: D{P, Q) = {<o{z{P) - z{Q))-^ + 0{l)]dz{P){dz{Q)f . From B{P, Q) 
and D[P^ Q) it is easy to construct meromorphic normalized (all a-periods vanish) 1-forms on C with 
single pole; in particular, if the pole coincides with ramification point P^, the natural local parameter 

Q=Pk 



near P^ is Xk{P) = ^/f[P) - A,,; then B{P,Pk) := g§ ^^^^ and D{P,Pu) :- js^^j^ 
meromorphic normalized 1-forms on L with single pole at P^ and the following singular parts: 

^(^' = + \sB{Pk) + 0(1)} dxk{P) ; D{P, Pk) = [j^^^^ + 0(1)} dxkiP) (7) 

as P ^ Pfc, where SBiPk) is the Bergmann projective connection computed at the branch point Pk 
with respect to the local parameter Xk{P)- 

Equations © in order look as follows (we write only equations with respect to Vi): 

r rpl 

-Y\P) (8) 



mf{p)) 



where the is the (taken with minus sign) l/N"^ contribution to the resolvent W. The function 
can be computed using the loop equations j3j and the "normalization conditions" 

i Y\P)df{P) = (9) 

over all basic a-cycles (these conditions mean that the "filling fractions" do not have the 1 /N"^ correc- 
tion). We introduce also the polynomial 

£{x,y) := {Vi{x) - y){V2{y) - x) - V{x,y) + 1 , (10) 

the function U{x, y), which is a polynomial in y and rational function in x: 

and rational function U{x, y, z): 

Then the loop equation looks as follows: 

„) = . - vAy) + - - ^^^^ . (13) 

y — Y [x) N'^y—Y{x) 
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where Y{x) := Vi'{x) — W(x); it arises as a corollary of reparametrization invariance of the matrix 
integral © 0. To use the loop equation effectively we need to consider the expansion of all of 

their ingredients. In this way we get the following expression for Y^: 



Y\P)df{P) 



V\f{P),g{P))df{P) 
£^Jf{P),g{P)) 



+ 



B{P,Q) 



Q^P ■■ f{Q)=f{P) 



df{Q) g{P)-g{Q) 



(14) 



where £y{x,y) means partial derivative with respect to the second argument. All ingredients of ()14() 
arise already in the leading term, except V^. However, from (jSJ we see that V{x,y) is a polynomial 
of degree di — 1 with respect to x and d2 — 1 with respect to y; moreover, the coefficient in front of 
^di-iyd2-i Joes not have l/N"^ correction. Thus we can conclude that the one-form Y^{P)df{P) is 
non-singular on the spectral curve outside of the branch points P^ (where it has poles of order 4); 
moreover, the first term in (|14() is non-singular on C (the first order zeros of £y at the branch points 
are cancelled by first order zeros of df[P) at these points). The form of singular parts at Pm allows 
to determine Y^(P)df{P) completely in terms of differentials B{P,Pf^) and D{P,Pf:) if we take into 
account the absence of correction to the "filling fractions" Q; the result looks as follows: 



mi 



y«(p)d/(p) = 



k=l 
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96g'iPk) 



D{P,Pk) + 



g"'{Pk) SBjPk] 
965'2(Pfc) 2Ag'{Pk 



B{P,Pk) 



(15) 



Then the solution of Q, (|15|) . which is symmetric with respect to the projection change (and, therefore, 
satisfies also equations ^ with respect to V2), looks as follows: 



if 1 ™i 1 

= -\nUf{v,,+,f--^'^dg{Pk) \ + 



d2 + 3 
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In do 



(16) 



where r/ is the so-called Bergmann tau- function on Hurwitz space [7j, which satisfies the following 
system of equations with respect to the branch points A^: 



Alnr, = -l5MP.); 



(17) 



In derivation of (|15|) we have used the following variational formulas, which can be easily proved in 
analogy to Rauch variational formulas: 



5\k 



mf{p)) 



g'{Pk)dfiP)=B{P,Pk) , 



(18) 



Hg'iPk)} 



mf{P)) 



dfiP) = - 
f(P) ' 4 



D{P,Pk) 



9"'{Pk) 



B{p,Pk: 



(19) 



g'iPk) 

The Bergmann tau-function (|17|) appears in many important problems: it coincides with isomon- 
odromic tau-function of Hurwitz Frobenius manifolds [SI, and gives the main contribution to G- 
function (solution of Getzler equation) of these Frobenius manifolds; it gives the most non-trivial term 
in isomonodromic tau-function of Riemann-Hilbert problem with quasi-permutation monodromies. Fi- 
nally, its modulus square essentially coincides with determinants of Laplace operator in metrics with 
conical singularities over Riemann surfaces [7j. The solution of the system ()17() looks as follows j7j- 
Define the divisor (df) = — 2oo/ — {d2 + l)oOg + YlT=i P^ ■= Yl^lf^ r^Dk- Choose some initial point 
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P ^ C and introduce corresponding vector of Riemann constants and Abel map Aa{Q) = /p Wa 
{wa form the basis of normalized holomorphic 1-forms on C). Since some points of divisor {df) have 
multiplicity 1, we can always choose the fundamental cell C of the universal covering of C in such a way 
that A{{df)) = —2K^ (for an arbitrary choice of fundamental domain these two vectors coincide only 
up to an integer combination of periods of holomorphic differentials), where the Abel map is computed 
along the path which does not intersect the boundary of C. The main ingredient of the Bergmann 
tau-function is the following holomorphic multivalued (1 — g)(7/2-differential C{P) (the higher genus 
analog of Dedekind eta- function) on C: 

where W{P) := deti<a,f3<g\\wp' ^\p)\\ denotes the Wronskian determinant of holomorphic differen- 
tials at point P; is the vector of Riemann constants with basepoint P. Introduce the quantity Q 
defined by 

m+N 

Q = [df{P)]'^C{P) n [E{P,D,)]^ ■ (21) 

k=l 

which is independent of the point P G C. Then the Bergmann tau-function p7p of Hurwitz space is 
given by the following expression [Jj: 

m+n 

Tf = Q'/^ n [E{Dk,Di)f-^ ; (22) 

k,l=l k<l 

together with (|16|) this gives the answer for l/N"^ correction in two-matrix model. If rj and Tg are 
Bergmann tau-functions H17() corresponding to divisors {df) and {dg), respectively, then 

" ^ „ (".i.+i)''^nt'y(Qt) , , 

where C = +^^^2'^^ 

is a constant independent of moduli parameters. Using the transformation 
(|2,'-{j) of the Bergmann tau-function under projection change, we find that the solution expression 
(|16() for satisfies also the necessary equations with respect to V2. This could be considered as a 
confirmation of consistency of our computation. Derivatives of function with respect to the 

filling fractions look as follows: 

^ = -i Y\P)df{P) ■ (24) 

these equations are counterparts of Seiberg-Witten type equations for (see for example 

One-matrix model. If potential V2 is quadratic, integration with respect to M2 in (pQ) can be 
performed explicitly, and the free energy ()16() gives rise to the free energy of one-matrix model. The 
spectral curve C in this case becomes hyperelliptic, and the formula (|16)) gives, using the expression 
for Tf obtained in j^l: 

F' = -ln{A^{detArllg'{X,)\ (25) 



A;=l 
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where X^, k = 1, . . . ,2g + 2 are branch points of A is their Wronskian determinant; A is the matrix 
of a-periods of non-normahzed holomorphic differentials on C; this agrees with previous results Hj. 

F^, isomonodromic tau-function and G -function of Frobenius manifolds. The genus 1 correction 
to free energy in topological field theories is given by so-called G-function, which for an arbitrary 
m-dimensional Frobenius manifold related to Hurwitz space looks as follows jHl IIUI Ej: exp{G} = 
Tj- ^^^rifcLil^sslpj. ^}~^/^^, where Tf is the Bergmann tau-function, ip is an "admissible" one-form on 
underlying Riemann surface. If, trying to build an analogy with our formula (|16() for F^, we formally 
choose ^p{P) = dg{P), the formula for G-function coincides with (|16)) up to small details like sign, 
additive constant, and the highest coefficient of polynomial ¥2- However, the differential dg is not 
admissible, and, therefore, does not really correspond to any Frobenius manifold; therefore, the true 
origin of his analogy is still unclear 

F^ and determinant of Laplace operator. Existence of close relationship between F^ and de- 
terminant of Laplace operator was suggested by several authors (see e.g. for hermitial one- 
matrix model, for hermitian two-matrix model and, finally, |14j for normal two-matrix model with 
simply-connected support of eigenvalues). After an appropriate regularization the (formal) determi- 
nant of Laplace operator A-^ over C in the singular metric \df{P)\'^ is given by the expression 
detA-^ = C^{det3=B}|T/|2 , where A is a regularized area of L, B is the matrix of 6-periods of L, C 
is a constant. In the "physical" case of real coefficients of Vi and V2 and real filling fractions the em- 
pirical expression for InjdetA-^} differs from our expression (|16|) by several explicit terms. Therefore, 
the relationship between hermitian and normal two-matrix models on the level of F^ seems to be not 
as straightforward as on the level of functions F^ • 
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